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FOREWORD  J 


This  final  report  describes  the  work  accomplished  by  Computer 
Sciences  Corporation  for  the  Aircraft  and  Crew  Systems  Technology 
Directorate  of  the  Naval  Air  Development  Center,  Warminster, 
Pennsylvania,  in  accordance  with  the  requirements  defined  by  Task 
Order  0047,  "escape  System  Trajectory  Sensitivity  Analysis", 
issued  under  Contract  1.62269-75-C-0001.  This  work  was  sponsored  by  the 
Naval  Air  Systems  Command  under  AirTask  A5365360/Q0 1D/3W0625-0000 
and  was  monitored  by  Mr.  Louis  D'Aulerio  and  Mr.  Alan  Cantor  of 
NADC. 
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- SECTION  Y 


At  the  request  of  the  Naval  Air  Development  Center  Crew  Systems  Depart- 
ment, the  Computer  Sciences  Corporation  initiated  a study  to  identify 
significant  input  parameters  to  two  existing  mathematical' 'seat  ejection 
models  (DEDALUS,  a two  dimensional  model  and  ICARUS,  a three  dimensional 
model) . CSC  was  also  requested  to  ascertain  the  models ' response  to 
variations  in  the  incut  parameter  values. 


1.1  INPUT  DATA  SELECTION 

Two  physically  different  seat  ejection  configurations  (rocket  position, 
rocket  angle,  seat/man  center  of  gravity,  etc.)  were  used  in  the  study. 
The  following  eight  input  parameters  were  selected  to  be.  varied:  the 
seat/man  center  of  gravity,  catapult  thrust,  rocket  thrust,  rocket 
angle,  rocket  ignition,  rocket  position,  seat/man  moments  of  inertia 
and  the  aerodynamic  coefficients.  The  above  eight  parameters  are  common 
to  both  models. 


1.2  METHODOLOGY  USED  IN  STUDY 

The  effects  on  the  seat/man  ejection  trajectory,  by  varying  the 
parameters  by  a plus  or  minus  twenty  percent  for  both  configurations, 
were  analyzed  under  zero  and  six  hundred  knots  initial  velocity  con- 
ditions. The  varied  parameter  ejection  trajectories  and  their  corre- 
sponding base  case  trajectory  were  then  plotted.  The  mean  values  and 
standard  deviations  of  the  differences  between  the  varied  parameter's 
trajectories  and  their  corresponding  base  case  trajectory,  over  spec- 
ified phases  of  the  trajectory  were  computed.  The  specified  phases 
were  catapult  ignition  to  catapult  separation,  rocket  ignition  to 
rocket  burnout,  rocket  burnout  to  drogue  chute  projection,  and  drogue 
projection  to  one  second.  The  parameters  were  then  ranked  using  their 
respective  means  and  standard  deviations  as  the  ranking  criterion.  The 
statistical  significance,  of  each  ranking  of  the  parameters  obtained, 
was  ascertained  by  computing  a Coefficient  of  Concordance,  the  associated 
F-test  value  at  the  95%  confidence  level  and  Spearman's  Rank  Correlation 
Coefficient.  These  coefficients  were  computed  as  a function  of  the  above 
rankings. 


1.3  ANALYSIS  OF  RESULTS  OBTAINED 

Analysis  of  the  results  obtained  from  the  study  indicates  that  all 
eight  variables  are  significant.  Each  of  the  eight  parameters,  when 
varied  by  a plus  or  minus  twenty  percent,  is  ranked  either  first,  second 
or  third  at  least  once  with  respect  to  the  amount  it  caused  a trajectory 
component  to  differ  from  the  component ' s corresponding  base  case  set 
of  values. 


The  following  factors  were  found  to  influence  the  individual  and  com- 
bined trajectory  component  rankings  obtained  during  the  study: 


1)  the  trajectory  components  responded  differently  from  one 
another  to  the  variations  in  the  input  parameter  values. 

2)  the  rankings  of  the  parameters  are  effected  by  initial 
velocity  differences  (zero  knots  versus  600  knots) . 

3)  the  rankings  of  the  parameters  are  effected  by  the 
differences  existing  between  the  two  sets  of  base  case 
values. 

A separate  analysis,  not  included  in  the  above  main  analysis,  identified 
the  weight  of  the  seat/man  combination  as  being  a significant  input 
parameter. 

A limited  analysis  of  the  effects  on  the  ejection  trajectory  when 
two  or  three  of  the  input  parameters  are  varied  simultaneously  by  a 
plus  or  minus  twenty  percent  indicates  certain  combinations  of  the  para- 
meters tend  to  minimize  the  differences  between  the  base  case 
trajectory  and  the  varied  parameter's  trajectory,  while  other  com- 
binations tended  to  maximize  the  differences. 

An  additional  case  was  constructed  whereby  the  eight  input  parameters 
were  varied  by  a small  amount  simultaneously.  The  difference  between 
the  trajectory  generated  with  the  eight  parameters  set  at  these  values 
and  their  corresponding  base  case  trajectory  values  was  almost  as 
great  as  the  maximum  difference  observed  when  the  eight  parameters 
were  varied  one  at  a time  by  a plus  or  minus  twenty  percent. 

1.4  RECOMMENDATION 

A reccmmenda tion  as  to  how  the  model  (s)  could  be  modified  to  arrive 
at  an  envelope  of  trajectories  surrounding  the  base  case  trajectory 
is  given.  This  envelope  of  trajectories  would  be  generated  stochasti- 
cally taking  into  account  the  degree  of  uncertainty  associated  with 
specific  or  all  of  the  input  parameters  and  the  effects  of  varying 
the  parameters  simultaneously. 
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SECTION  2 

INTRODUCTION 


The  following  report  documents  the  work  completed  by  Computer  Sciences 
Corporation  under  Task  Order  Number  47,  Contract  N62269-75rC-0001  for 
the  Naval  Air  Development  Center,  Crew  Systems  Engineering  Division, 
Warminster,  Pennsylvania.  The  report  is  in  two  parts. 


2.1  DESCRIPTION  OF  THE  FIRST  PART  OF  THE  REPORT 

The  first  part  of  the  report,  Sections  3 to  6,  documents  the  methodology 
used  and  the  results  obtained  while  using  this  methodology  to  identify 
the  significant  input  parameters  to  two  existing  seat  ejection  mathema- 
tical simulation  models,  programmed  in  Fortran  and  executed  on  the  CDC 
6600  Computer. 


2.1.1  Description  of  the  Two  Seat  Ejection  Simulation  Models.  The  two 
models  referenced  here  are  DEDALUS,  a two  dimensional  model,  and  ICARUS, 
a three  dimensional  model.  DEDALUS  describes  the  seat/man  ejection  tra- 
jectory in  the  plane  of  the  horizontal  x-axis  and  vertical  z-axis.  The 
above  seat/man  spatial  coordinators  (x,  z)  are  computed  by  the  model 
at  each  integration  time  step.  The  angular  rate  of  rotation  (pitching 
rate)  in  this  plane  is  also  computed  by  the  model. 

The  three  dimensional  model  (ICARUS,  a six  degrees  of  freedom  model)  com- 
pletely describes  the  ejection  trajectory.  The  three  linear  spatial  coor- 
dinates (x,  y and  z)  and  the  three  angular  rates  of  motion  (p,  q and  r 
corresponding  to  roll,  pitch  and  yaw  respectively)  of  the  seat/man  ejec- 
tion trajectory  are  computed  by  the  model. 


2.1.2  Description  of  the  Base  Case  and  Varied  Parameter  Trajectories. 
As  an  aid  in  the  identification  of  the  significant  input  parameters  to 
the  ejection  seat  models,  it  was  decided  to  set  the  input  parameters 
to  their  base  case  values  and  to  generate  a seat/man  ejection  trajectory 
(via  the  models)  using  these  values  as  input  to  the  model  (s)  . The  base 
case  value  being  defined  here  to  mean  the  best  estimate  of  the  input 
parameter's  true  value  for  a given  physical  seat-ejection  system.  The 
input  parameters  were  varied  one  at  a time  and  an  ejection  trajectory  was 
generated,  via  the  model,  for  each  variable  varied. 


From  past  experience  in  the  use  of  the  models,  it  was  noted  that  when  an 
individual  parameter  was  varied  frcm  a given  initial  value,  the  variable 
had  to  be  varied  anywhere  from  fifteen  to  twenty  five  percent  before  a 
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consistent  noticeable  change  in  the  trajectory  (relative  to  the  initial 
trajectory)  was  observed.  For  the  purposes  of  this  study,  it  was  decided 
to  vary  the  input  parameters  one  at  a time  - first  by  a plus  twenty  per- 
cent then  by  a minus  twenty  percent.  An  ejection  trajectory  was  then 
generated  with  one  of  the  input  parameters  set  to  a plus  twenty  percent 
increase  in  value  over  its  base  case  value  while  the  other  input  para- 
meters were  held  constant  at  their  respective  base  case  values. 

Similarly,  an  ejection  trajectory  was  generated,  with  the  above  parameter 
set  to  a minus  twenty  percent  decrease  in  value  from  its  base  case  value. 
These  trajectories  were  then  used  to  identify  the  significant  input  para- 
meters and  to  measure  the  models'  response  (relative  to  the  base  case 
trajectory)  to  the  variation  in  the  parameter's  base  case  value. 


2.1.3  Sunmary  Report  Implementation.  In  order  to  reduce  the  time  and 
cost  of  generating  an  individual  ejection  trajectory  via  the  models  on 
the  computer,  a sunmary  report  displaying  the  above  linear  and  angular 
components  of  the  trajectory  was  designed  and  implemented  in  both  models. 
The  models  were  then  further  modified  to  generate  an  output  file  of  re- 
cords containing  the  time  followed  by  the  values  of  the  spatial  coordinates 
and  the  angular  rate(s)  of  the  trajectory  corresponding  to  the  time 
value  indicated.  Thus,  a separate  file  containing  the  profile  cf  time  and 
the  trajectory  component  values  was  retained  for  each  trajectory  generated 
by  the  model (s) . 


2.1.4  Description  of  the  Plotting  and  Statistical  Programs  used  in  the 
Study.  The  above  files  were  later  used  as  input  to  two  other  programs. 
One  program,  based  upon  the  information  contained  in  the  files  caused 
the  trajectory  data  to  be  plotted.  The  other  program,  which 
was  especially  written  for 'this  study,  performed  a statistical  comparison 
between  any  two  of  the  above  generated  trajectory  files.  This  program 
was  used  to  perform  a statistical  comparison  between  a given  base  case 
trajectory  component  profile  and  a corresponding  varied  parameter  tra- 
jectory component  profile.  A detailed  explanation  of  how  the  above 
generated  trajectories  were  used  to  identify  the  significant  input 
parameters  and  to  measure  quantitatively  the  models'  response  to  an 
input  parameter  being  varied  is  contained  in  the  succeeding  sections 
of  this  report. 


2.1.5  Explanation  of  the  Numerical  Integrating  Algorithm  Modification. 
Also,  to  further  reduce  the  time  required  and  the  associated  cost  of 
generating  an  ejection  trajectory  using  the  three  dimensional  model 
(ICARUS) , the  fourth  order  Runge-Kutta  integrating  algorithm  used  in 
ICARUS  was  replaced  with  a Runge-Kutta/Simpson ' s rule  integrating  scheme. 
The  Runge-Kutta  algorithm  was  used  to  obtain  starting  solutions  for  the 
Simpson' s rule  algorithm  and  to  ascertain  the  precise  time  a position 
or  time  event  occurred. 
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The  Simpson's  rule  algorithm  was  used  for  all  constant  step  size  inte- 
grations once  the  required  number  of  starting  solutions  were  obtained. 
The  replacement  of  the  Runge-Kutta  algorithm  with  the  Runge-Kutta/ 
Simpson's  Rule  scheme  resulted  in  reducing  the  cost  of  generating  a 
seat  ejection  trajectory  on  the  computer  by  a factor  of  one-half. 

A comparison  of  the  trajectory  components  generated  using  the  Runge- 
Kutta  throughout  the  trajectory,  as  compared  to  the  Runge-Kutta/Simpson 
rule  scheme,  indicated  that  for  the  low  velocity  cases  - the  trajectory 
components  generated  by  both  integrating  schemes  were  within  two  per- 
cent of  each  other  after  four  seconds  of  the  trajectory  simulation. 
However,  a comparison  of  the  trajectory  components  generated  by  both 
schemes  at  the  600  knot  initial  velocity  indicated  that,  after  1.5 
seconds  of  the  trajectory,  the  components  disagreed  by  as  much  as  five 
percent.  In  addition,  both  integrating  schemes  exhibited  instability 
when  a integration  step  size  of  .005  or  less  was  used  after  drogue 
chute  projection. 

Based  upon  the  above  integrating  cost  analysis  and  step  size  require- 
ments, it  was  decided  to  use  the  Runge-Kutta/Simpson  rule  algorithm, 
keeping  the  integrating  step  size  constant  at  .001  throughout  the 
trajectory  simulation,  and  to  perform  the  input  parameter  study  using 
the  values  obtained  for  one  second  of  the  simulation  trajectory;  thus, 
minimizing  the  errors  due  to  the  numerical  integration. 


2.1.6  Preliminary  Case  Studies.  Finally,  a series  of  preliminary  case 
studies  were  executed  on  the  computer  via  the  two  models.  Because  of 
the  rather  low  cost  per  an  ejection  trajectory  generation  associated 
with  the  two  dimensional  model,  DEDALUS  was  used  mainly  during  this 
initial  stage  of  the  study. 

The  purpose  of  this  preliminary  investigation  was  to: 

1)  insure  that  a + 20%  variation  level  was  sufficient  to 
detect  differences  between  the  base  case  and  the  varied 
parameter  trajectory  components. 

2)  insure  that  both  models  responded  similarly  to  the  var- 
iables varied  taking  into  account  one  was  a two  dimensional 
and  the  other  a three  dimensional  model. 

3)  perform  a preliminary  screening  as  to  which  of  the  input 
parameters  were  to  be  used  in  the  main  body  of  the  study. 
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2.1.7  Results  of  the  Preliininary  Case  Studies . As  a result  of  the 
above  preliminary  study,  the  following  was  verified  or  accomplished: 

1)  a + 20%  variation  level  was  sufficient  to  detect  differ- 
ences between  the  two  trajectories. 

2)  both  models  do  respond  in  the  same  way  to  the  variables 
varied  taking  into  account  one  is  a 2 -dimensional  model 
while  the  other  is  a 3-dimensional  model. 


3)  a preliminary  screening  of  the  input  parameters  yielded 
the  following  list  of  eight  parameters  to  be  used  in  the 
main  portion  of  the  study: 

a)  the  seat/man  center  of  gravity  (S/M  C.G.) 

b)  the  catapult  time/thrust  table  (Cat  Thr) 

c)  the  rocket  time/thrust  table  (Roc  Thr) 

d)  the  distance  in  the  vertical  direction  that  the 
rocket  must  travel  (relative  to  its  initial  pos- 
ition) before  it  is  ignited  (Roc  On) 

e)  the  seat/man  moments  of  inertia  (S/M  M.I.) 

f)  the  rocket  location  with  respect  to  the  catapult 
(Roc  Pos) 

g)  the  seat/man  aerodynamic  coefficient  table  (Aero) 

Section  3 of  this  report  gives  a detailed  description  of  the  above  eight 
parameters. 


Section  4 explains  the  methodology  used  to  quantify  and  to  generate  the 
rankings  of  the  eight  input  parameters  with  respect  to  the  amount  they 
disturbed  the  six  components  of  the  trajectory  (relative  to  their  re- 
spective base  case  values)  when  varied. 

Section  5 explains  the  results  obtained  from  the  study  while  using  the 
above  methodology. 

Section  6 recommends  an  approach  to  be  taken  in  determining  the  model  (s) 
response  to  the  simultaneous  varying  of  certain  or  all  of  the  input 
parameters  from  their  base  case  values. 


2.2  DESCRIPTION  OF  THE  SECOND  PART  OF  THE  REPORT 

The  second  part  of  the  report.  Section  7,  documents  the  work  completed 
to  date  by  Computer  Sciences  Corporation  in  modifying  and  extending  the 
capabilities  of  the  above  two  seat  ejection  models . 
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SECTION  3 j 
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DEFINITION  OF  THE  EIGOT  INPOT  PARAMETERS 
SELECTED  TO  BE  VARIED  IN  THE  STUDY 


3.1  SEAT/MAN  CENTER  OF  GRAVITY  (S/M  C.G.)  . 

Center  of  gravity  of  the  seat/man  combination.  The  x,  y and  z 
coordinates  of  this  parameter  being  specified  in  the  Catapult  Coordi- 
nate System  (CCS)  in  units  of  feet. 


3.2  CATAPULT  TIME/THRUST  TABLE  (CAT  THR) . 

Table  of  values  displaying  the  elapsed  time  relative  to  catapult  igni- 
tion and  the  corresponding  force  that  the  catapult  exerts  on  the  seat/ 
man  combination  at  that  specific  time  in  units  of  lbs. 


3.3  ROCKET  TIME/THRUST  TABLE  (ROC  THR)  . 

Table  of  values  displaying  the  elapsed  time  relative  to  rocket  ignition 
and  the  corresponding  force  that  the  rocket  exerts  on  the  seat/man 
combination  at  that  specific  time  in  units  of  lbs. 


3.4  ROCKET  THRUST  ANGLES  (ROC  ANG) . 

The  angles  that  the  rocket  thrust  line  makes  with  the  catapult  coordi- 
nate axis  in  units  of  degrees. 


3.5  ROCKET  IGNITION  LENGTH  (ROC  ON). 

The  distance  the  rocket  must  travel  in  the  vertical  direction  relative 
to  its  static  position  before  the  rocket  will  ignite.  The  distance 
being  measured  in  units  of  feet. 


3.6  ROCKET  POSITION  (ROC  POS) . 

Position  of  the  rocket  relative  to  the  catapult.  The  x,  y and  z 
coordinates  being  given  in  the  Catapult  Coordinate  System  (CCS)  in 
units  of  feet. 


3.7  SEAT/MAN  MOMENTS  OF  INERTIA  (S/M  M.I.). 

The  three  principal  ocmponents  of  the  seat/man  moments  of  inertia  (Ixx, 
Iyy,  Izz)  in  units  of  slug-feet**2. 
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3.8  AERODYNAMIC  COEFFICIENTS  (AERO). 

The  aerodynamic  force  and  moment  coefficients  used  to  arrive  at  the 
forces  and  moments  acting  on  the  seat/man  combination.  Coefficients 
derived  as  a function  of  the  angle  of  attack,  Mach  number  and  sideslip 
(3-dimensional  model  only) . 
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SECTION  4 

EXPLANATION  OF  THE  METHODOIOGY  USED  IN  STUDY 


The  following  describes  the  methodology  used  and  results  obtained  while 
using  this  methodology  to  identify  the  significant  input  parameters  to 
two  existing  seat  ejection  mathematical  simulation  models.  The  method- 
ology was  also  used  to  quantitatively  measure  the  models'  response  to 
variations  in  selected  input  parameter  values  frcm  their  initial  (base 
case)  values. 


4.1  GENERATION  OF  THE  BASE  CASE  AND  VARIED  PARAMETER  TRAJECTORIES 

Two  physically  different  sets  of  input  parameters  (rocket  position, 
seat/man  center  of  gravity,  moments  of  inertia,  etc.)  were  used  in  the 
study  and  are  referenced  in  this  report  as  Case  1 and  Case  2. 

With  the  input  parameters  set  to  their  corresponding  base  case  values,  two 
ejection  trajectories  were  generated  for  each  of  the  above  cases  by 
the  model (s).  One  ejection  trajectory  was  generated  simulating  a zero 
altitude-zero  knot  initial  velocity  ejection,  the  other  simulating  a 
zero  altitude-600  knots  initial  velocity  ejection.  In  this  study,  these 
four  trajectories  are  referenced  as  the  base  case  trajectories;  namely, 
Case  1 at  zero  knot  initial  velocity,  Case  1 at  600-knots,  Case  2 at 
zero  knot  and  Case  2 at  600-knots. 

The  eight  input  parameters,  selected  to  be  varied  in  this  study,  were 
then  varied  individually.  Each  variable  was  first  varied  from  its 
base  case  value  by  a plus  twenty  percent  then  by  a minus  twenty  percent, 
while  the  other  input  parameters  were  held  constant  at  their  initial 
base  case  values.  An  ejection  trajectory  was  generated,  via  the  3- 
dimensional  model  ICARUS,  for  each  case,  velocity,  percent  variation 
and  varied  parameter  combination  - resulting  in  64  ejection  trajectories 
being  generated  (2X2X2X8) . 


4.2  FORMAT  OF  THE  FILES  CONTAINING  THE  TRAJECTORIES 

The  components  of  the  above  64  varied  parameter  trajectories  plus  the  4 
base  case  trajectories  generated  by  the  model  were  retained  on  seoarate 
files  in  the  computer . The  record  format  of  the  files  created  being: 
time,  followed  by  the  3-linear  spatial  coordinate  (x,  y and  z)  values 
of  the  seat/man  center  of  gravity  with  respect  to  an  earth  fixed  coordi- 
nate system  and  these  followed  by  the  3-angular  rates  of  rotation  (p,  q 
and  r)  of  the  seat/man  combination  about  a seat/man  coordinate  system 
with  origin  at  the  seat/ran  center  of  gravity.  Each  file  containing 
approximately  one  thousand  records,  one  record  for  each  time  step 
(.001  secs.)  of  the  simulated  trajectory. 
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4.3  USE  OF  THE  FITES  TO  GENERATE  PLOTS  AND  THE  RESIDUAL  VALUES  OF 

THE  TRAJECTORY  COMPONENTS 

The  above  files  were  then  used  as  input  data  to  two  separate  Drograms, 
one  program  being  a plot  program.  This  program  was  used  to  generate 
plots  of  the  vertical  and  horizontal  coordinates  of  the  seat/man 
center  of  gravity.  Appendix  A contains  the  plots  of  the  Case  1 varied 
parameter  trajectories  versus  their  respective  base  case  trajectory. 

The  files  were  also  processed  by  a second  program  which  was  especially 
written  for  this  study.  This  program  was  used  to  compute  the  difference 
(residual  value)  between  a given  base  case  trajectory  component  profile 
and  a corresponding  varied  parameter  profile.  The  six  components  (x,  y, 
z,  p,  q and  r)  of  the  base  case  trajectory  profile  being  compared  indi- 
vidually to  their  corresponding  components  contained  in  the  varied  para- 
meter profile.  The  difference  between  the  six  component  pairs  being 
computed  over  the  1 second  time  profile. 

(xbase,i  -xvar  par,i,  ^base , i "/var  par,i,  ...,  rbase,i 
-rvar  par,i  where  i = 0,  1,  2,  ...,  1000  corresponding 
to  the  trajectory  simulation  time  starting  at  zero,  catapult 
ignition,  being  incremented  every  .001  secs,  and  terminating 
1 sec.  later) . 

The  mean  values  and  standard  deviations  of  the  six  residual  value  vectors 
were  then  computed  over  selected  intervals  of  the  vectors.  The  intervals 
corresponding  to  the  following  phases  of  the  seat/man  ejection  simulation: 

1)  Catapult  ignition  (time  zero)  to  catapult  separation. 

2)  Rocket  ignition  to  rocket  burnout. 

3)  Rocket  burnout  to  drogue  chute  projection. 

4)  Drogue  chute  projection  to  1 second  of  the  simulation. 

Using  the  corresponding  mean  value,  standard  deviation  and  the  values 
of  the  residuals  within  a given  phase  of  the  trajectory,  the  residual 
value  for  each  phase,  trajectory  component  combination  was  then  com- 
puted at  the  95%  confidence  level  |Rgc|  = |R  + a p | . The  value  of  the 
residual  at  the  95%  confidence  level  being  defined9nere  to  mean  that 
value  of  the  residual  which  is  95%  of  the  time  greater  than  or  equal 
to  the  value  of  the  residuals  contained  within  the  phase  or  conversely 
is  exceeded  in  value  5%  of  the  time  by  the  residual  values  contained 
in  the  phase. 

The  residual  values  at  the  95%  confidence  level  oomputed  above  were 
then  tabulated  by  case,  velocity,  percent  variation,  trajectory  component 
and  phase.  Appendix  B displays  the  tabulated  values  of  the  residuals  at 
the  95%  level.  Appendix  E presents  in  mathematical  form  the  equations 
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used  to  compute  the  above  residual  values. 


4.4  USE  OF  TOE  RESIDUAL  VALUES  TO  RANK  THE  EIGHT  INPUT  PARAMETERS 

The  tabulated  residual  values  were  then  used  to  rank  the  eight  input 
parameters  with  respect  to  the  amount  each  parameter  when  varied  by 
a plus  or  minus  twenty  percent  caused  the  seat/man  ejection  trajectory 
to  differ  from  its  corresponding  base  case  trajectory.  An  individual 
ranking  of  the  eight  parameters  was  determined  for  each  case,  velocity, 
percent  variation,  trajectory  component  and  phase  combination. 

The  procedure  used  to  form  the  individual  rankings  of  the  eight  para- 
meters was  as  follows: 

1)  For  a specific  case,  velocity,  percent  variation,  trajectory 
component  and  phase  combination  - the  parameter  associated 
with  the  largest  residual  value  - was  ranked  first  and 
assigned  the  numerical  value  of  1.  Then  the  parameter 
associated  with  the  second  largest  residual  value  was 
ranked  second  and  assigned  the  numerical  value  of  2.  This 
procedure  was  continued  until  all  eight  parameters  were 
ranked  and  assigned  a numerical  value  corresponding  to 
their  relative  position  in  the  list. 

2)  This  procedure  was  repeated  for  every  case  (2) , velocity 
(2),  percent  variation  (2),  trajectory  component  (6)  and 
phase  combination  (3)  yielding  108  (2X2X2X6X3)  numerical 
rankings  of  the  eight  input  parameters. 

The  catapult  phase  of  the  trajectory  was  excluded  from  the  above  rank- 
ings procedure  as  the  rocket  input  parameters  are  not  used  in  the  com- 
putation of  trajectory  until  the  very  end  of  the  phase. 


4.5  COMPUTATION  OF  THE  COEFFICIENT  OF  CONCORDANCE  AND  SPEARMAN'S 

RANK  CORRELATION  COEFFICIENT 

The  following  correlation  coefficients  were  then  computed  based  upon 
these  numerical  rankings  of  the  eight  input  parameters. 

A Coefficient  of  Concordance  was  computed  to  measure  the  degree  of 
overall  agreement  existing  between  the  phase  rankings  of  the  variables 
for  every  case,  velocity,  percent  variation  and  trajectory  component 
combination.  The  possible  values  for  the  coefficient  computed  being 
in  the  interval  of  zero  to  one.  A zero  value  obtained  for  the  Concor- 
dance Coefficient  indicating  that  the  overall  agreement  existing  among 
the  phases  as  to  the  ranking  of  the  parameters  is  no  better  than  a 
random  ranking  of  the  parameters.  Whereas,  a value  of  one  computed 
for  the  coefficient  indicates  that  the  phase  rankings  are  in  complete 
agreement  with  one  another. 
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Snedecor's  F-test  both  at  the  95%  and  90%  level  was  used  to  test  the 
statistical  significance  of  the  Concordance  Coefficients  computed. 

That  is,  an  F-value  to  be  tested  was  computed  for  each  of  the  above 
Coefficients  of  Concordance  derived  from  the  numerical  rankings. 

The  F-value  being  computed  as  a function  of  the  number  of  trajectory 
phase  rankings  combined  (pooled)  and  their  corresponding  computed 
Coefficient  of  Concordance.  The  degrees  of  freedom  to  be  used  in 
entering  Snedecor's  F-tables  were  then  computed  as  a function  of 
the  number  of  trajectory  phase  rankings  combined  and  the  number  of 
items  to  be  ranked. 

Entering  Snedecor's  tables  of  F with  the  above  derived  degrees  of 
freedom  for  the  greater  and  lesser  ratio  estimate,  each  of  the  above 
computed  F-values  were  tested  at  both  the  90%  and  95%  confidence  levels. 

If  the  value  of  F found  in  the  tables  was  less  than  the  F-value  com- 
puted - than  the  interpretation  of  the  Coefficient  fo  Concordance  was 
significant  at  the  confidence  level  indicated  for  the  table  of  values 
entered. 

At  the  95%  level,  the  interpretation  of  the  coefficient  being  rather 
significant  statistically.  That  is,  the  Concordance  Coefficient  being 
tested  - does  truly  indicate  the  degree  of  correlation  existing  among 
the  phases  in  ranking  the  parameters, 

Spearman's  Rank  Correlation  Coefficient  was  then  computed  to  measure 
the  degree  of  agreement  existing  between  any  two  of  the  above  phases 
of  the  trajectory.  The  range  of  possible  values  for  this  coefficient  being 
a plus  or  minus  one.  A minus  one  indicating  that  between  any  two 
specific  phase  rankings  of  the  parameters  - the  degree  of  agreement  is 
negatively  correlated  - that  is,  the  rankings  of  the  parameters  are  completely 
the  opposite  of  one  another.  A zero  value  indicates  there  is  no  correlation 
in  the  rankings  between  the  two  phases  and  finally  a positive  one 
indicating  complete  agreement.  This  coefficient  was  then  used  to 
determine  the  amount  of  agreement  existing  between  the  phases  in 
forming  a new  ranking  of  the  parameters  based  upon  the  combined  numer- 
ical values  of  the  individual  phase  rankings  of  the  parameters.  The 
individual  component  rankings  of  the  parameters  were  then  pooled  together 
(within  case  and  between  cases)  and  new  rankings  of  the  variables  were 
formed  based  upon  these  poolings.  The  corresponding  Coefficient  of 
Concordance,  F-values  and  Rank  Correlation  Coefficients  were  computed 
for  each  pooling  combination. 

Appendix  C presents  the  tabulated  values  of  the  above  coefficients  and 
F-values  derived.  The  individual  and  pooled  rankings  of  the  eight 
parameters  is  also  displayed  in  this  appendix.  Appendix  F contains 
a mathematical  presentation  of  the  ranking  procedure  used  to  derive 
the  Coefficient  of  Concordance  and  Spearman's  Rank  Correlation  Coefficients. 
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SECTION  5 

ANALYSIS  OF  THE  COEFFICIENTS  AND  RANKINGS  OF  THE 
INPOT  PARAMETERS  OBTAINED  FROM  THE  STUDY 


Appendix  C (Tables  1 thru  9,  Rl  thru  R12)  present  the  tabulated  values  of 
the  Coefficients  of  Concordance,  their  corresponding  F-values,  Spearman's 
Rank  Correlation  Coefficients  and  the  various  rankings  of  the  eight  input 
parameters  obtained  during  each  step  of  the  study.  As  presented  in 
Section  4 of  this  report,  these  coefficients  were  computed  as  a function 
of  the  (individual  and  pooled)  trajectory  component  numerical  rankings 
of  the  parameters. 


The  numerical  rankings  of  the  parameters  were  ascertained  by  the  amount 
each  input  parameter  (when  varied)  disturbed  the  values  of  the  trajectory 
components  from  their  corresponding  base  case  values. 


5.1  INDIVIDUAL  TRAJECTORY  COMPONENT  RANKINGS  AND  COEFFICIENTS 

Tables  1 to  4 present  the  tabulation  (by  case,  velocity  and  percent  var- 
iation) of  the  correlation  coefficients  obtained  from  the  individual 
trajectory  component  rankings  of  the  eight  input  parameters.  Tables 
Rl  to  R4  display  the  actual  ordering  of  the  parameters  obtained  from  the 
individual  numerical  rankings. 

An  examination  of  the  values  contained  in  the  Tables  1 thru  4 indicates 
there  exists  a relatively  high  correlation  among  the  trajectory  phases 
in  the  ranking  of  the  parameters  by  the  individual  trajectory  components. 

In  addition,  there  is  a 95%  probability,  except  for  two  instances,  that 
the  overall  degree  of  agreement  indicated  by  the  Concordance  Coefficients 
can  be  accepted  as  being  valid.  The  two  exceptions  being: 

1)  the  numerical  rankings  assigned  to  the  input  parameters 

by  the  z coordinate  of  the  trajectory  in  the  -20%  variation 
section  of  Table  2 resulted  in  a computed  F-value  of  2.55. 

This  value  was  less  than  the  F-table  value  of  2.92  required  for 
the  acceptance  of  the  Coefficient  of  Concordance  value  (.56)  at 
the  95%  confidence  level.  However,  the  coefficient  can  be  accepted 
at  the  90%  confidence  level.  The  F-table  value  at  this  level 
being  2.28. 

2)  the  numerical  rankings  assigned  to  the  input  parameters 
by  the  y coordinate  of  the  trajectory  in  the  +20%  varia- 
tion section  of  Table  2 resulted  in  a computed  F-value 

of  1.44.  This  value  of  F was  found  not  to  be  significant 
at  either  the  95%  or  90%  confidence  level. 

Finally,  the  relatively  high  values  of  the  Spearman's  Rank  Correlation 
Coefficients,  presented  in  the  Tables  1 thru  4,  indicates  that  the  new 
rankings  formed  by  combining  the  individual  phase  rankings  of  the  para- 
meters are  statistically  significant. 
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5.2  WITHIN  CASE  POOLED  RANKINGS 

Tables  5 and  6 present  the  correlation  coefficients  and  corresponding 
F-values  by  case,  velocity  and  percent  variation.  The  values  displayed 
in  the  tables  were  obtained  by  the  pooling  of  the  individual  linear 
(x,  y and  z combined)  and  angular  rate  (p,  q and  r combined)  component 
rankings  of  the  input  parameters.  The  upper  and  middle  sections  of 
Table  5 (Case  1)  and  Table  6 (Case  2)  present  the  coefficients  obtained 
by  the  within  Case  1 and  within  Case  2 pooling  of  the  linear  and  angular 
rankings . The  lower  section  (0  +600)  of  the  tables  present  the  coef- 
ficients obtained  by  the  further  within  case  combining  of  the  linear 
and  angular  rankings,  this  time  across  velocities  (zero  knot  combined 
with  600  knots  initial  velocity) . 

A ccnparison  between  the  values  displayed  in  the  upper  and  middle  sec- 
tions of  Tables  5 and  6 and  their  corresponding  individual  component 
values  contained  in  Tables  1 thru  4 indicates,  as  the  component  rankings 
are  grouped  together,  the  amount  of  agreement  as  to  the  ranking  of  the 
eight  input  variables  tends  to  decrease.  That  is,  the  coefficients  of 
Concordance  and  the  Rank  Correlation  Coefficients  derived  from  the 
within  care,  within  velocity  pooling  of  the  linear  and  angular  rankings 
tend  to  be  lower  in  value  than  their  corresponding  individual  component 
coefficient  values  found  in  Tables  1 thru  4.  This  decrease  in  the  values 
of  the  coefficients,  when  the  trajectory  component  rankings  are  pooled, 
indicates  that  the  individual  component  rankings  of  the  eight  input 
parameters  are  slightly  different  from  one  another. 

The  lower  section  of  Tables  5 and  6 present  the  values  of  the  Coefficients 
obtained  when  the  linear  and  angular  rankings  were  further  pooled  across 
velocities  (zero  velocity  combined  with  600  knots  velocity) . A comparison 
between  the  values  presented  in  the  lower  section  of  Tables  5 and  6 to 
their  corresponding  values  contained  in  the  upper  and  middle  sections  of 
the  same  tables  indicates,  the  linear  rankings  when  further  pooled  across 
velocities  results  in  the  degree  of  agreement  in  the  ranking  of  the 
parameters  being  lowered . This  was  true  for  both  Case  1 and  Case  2 coolings. 
However,  a comparison  of  the  angular  poolings  across  velocities  to  their 
corresponding  component  poolings  indicates  the  following: 

1)  the  combining  of  the  Case  l's  3 angular  rankings  at 
zero  velocity  with  Case  l's  3-angular  rankings  at  600 
knots  results  in  the  further  reduction  of  the  correlation 
coefficient  values  relative  to  the  values  obtained  for  the 
two  corresponding  within  velocity  poolings  of  the  coef- 
ficients. 

2)  the  combining  of  the  Case  2's  angular  rankings  across 
velocities  however,  results  in  obtaining  somewhat  of 
an  average  value  for  the  correlation  coefficients. 
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The  anomaly  between  the  Case  1 angular  component  pooling  across 
velocities  and  the  Case  2 angular  pooling  across  velocities  can  best 
be  explained  by  the  following  rational.  The  weight  of  the  seat/man 
combination  used  for  the  Case  1 at  600  knots  initial  velocity  was  20% 
lighter  than  the  weight  of  the  seat/ran  combination  used  for  Case  1 at 
zero  velocity  initial  condition.  A separate  sidebar  study  identified  the 
weight  of  the  seat/man  as  being  a significant  input  parameter.  Therefore, 
in  the  Case  1 pooling  of  the  angular  rates  across  velocities,  the  correla- 
tion coefficients  obtained  reflect  the  effect  of  the  difference  in  the 
seat/man  weight  along  with  the  difference  in  velocity  has  upon  the  pooled 
ranking  obtained.  The  weight  of  the  seat/man  was  the  same  in  Case  2 
for  both  the  zero  and  600  knots  initial  velocity  studies.  Thus,  the 
pooling  of  the  linear  rankings  accross  velocities  tends  to  indicate,  even 
with  the  difference  in  the  weight  of  the  seat/man  combination  effecting  the 
results  of  the  Case  1 poolings,  that  at  different  velocities  the  ranking 
of  the  input  parameters  can  be  different  frcm  on  another.  This  is  confirmed 
by  the  reduction  in  the  Case  2 pooled  velocity  correlation  coefficients 
obtained  in  this  case,  as  was  stated,  the  weights  were  the  same.  The 
values  of  the  Concordance  Coefficients  shown  in  Tables  5 and  6 were  all 
statistically  significant  at  the  90%  and  95%  confidence  levels. 

The  actual  rankings  of  the  eight  input  parameters  formed  at  each  step 
of  the  above  pooling  process  are  displayed  in  Tables  R5  thru  R8.  Multiple 
entries  in  a given  row/column  designated  cell  of  the  tables  indicates 
that  the  ranking  method  used  assigned  the  same  numerical  ranking  to 
the  parameters  displayed  in  the  cell.  As  stated  above,  the  pooling  of 
the  linear  and  angular  rankings  resulted  in  a reduction  in  the  degree 
of  agreement  in  the  new  ranking  of  the  parameters  formed,  and  when  the  rank- 
ings were  further  pooled  (zero  initial  velocity  combined  with  600  knots 
initial  velocity)  the  degree  of  agreement  was  reduced  further.  This 
suggests  the  actual  rankings  of  the  parameter  in  the  upper  and  lower  sec- 
tions of  Tables  F5  thru  R8  ( the  pooled  linear  and  pooled  angular  rankings) 
are  statistically  less  signif leant  than  the  rankings  of  the  parameters 
displayed  in  Tables  Rl  thru  R4  (the  individual  trajectory  component 
rankings) . Because  of  low  values  obtained  for  the  above  pooled  velocity 
correlation  coefficients,  the  corresponding  actual  rankings  of  the  para- 
meters are  not  displayed  in  the  tables. 


5,3  BETV2EEN  CASE  POOLED  RANKINGS 

Tables  7,  8 and  9,  R9  thru  R12  display  the  values  of  the  correlation 
coefficients  and  the  actual  rankings  of  the  eight  input  parameters 
obtained  by  pooling  Case  1 and  Case  2 individual  component  rankings. 
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5.3.1  Between  Case  Individual  Trajectory  Component  Rankings  of  the 
Parameters  Pooled.  Tables  7 and  8 present  the  values  of  the  coefficients 
obtained  by  combining  (within  velocity)  the  individual  trajectory  com- 
ponent rankings  of  Case  1 with  their  corresponding  component  ranidngs 
in  Case  2: 

(the  x-trajectory  component  rankings  of  Case  1 combined 
with  the  x-trajectory  component  rankings  of  Case  2,  the 
y-trajecotry  component  rankings  of  Case  1 combined  with 
the  y- trajectory  component  rankings  of  Case  2,  etc.) 

A comparison  of  the  coefficients  displayed  in  Table  7,  zero  velocity, 
to  their  corresponding  individual  component  values  contained  in  Tables 
1 and  2,  indicates  that  the  pooling  of  the  individual  linear  rankings 
(Case  1 components  combined  with  Case  2 components)  reduces  the  corre- 
lation coefficient  values  relative  to  their  corresponding  component 
values  found  in  Tables  1 and  2.  However,  a comparison  of  the  indivi- 
dually pooled  angular  rate  components  to  their  corresponding  components 
also  found  in  Tables  1 and  2 indicates  that  the  angular  coefficient 
remain  relatively  unchanged  except  for  the  effects  of  data  averaging 
which  is  inherent  in  the  ranking  procedure  used. 

A comparison  of  the  coefficients  displayed  in  Table  8 (obtained  by 
pooling  Case  1 and  Case  2 individual  linear  and  angular  rankings  of 
the  parameters  at  600  knots  initial  velocity)  to  their  corresponding 
individual  component  values  found  in  Tables  3 and  4 indicates  that 
both  the  linear  and  angular  coefficient  values  have  been  reduced  in 
value  from  their  corresponding  individual  component  values  found  in 
Tables  3 and  4. 

Thus,  at  the  low  velocity  initial  condition  there  is  evidence  that 
the  two  physically  different  sets  of  input  parameters  can  rank  the 
parameters  differently  as  evidenced  in  the  reduction  of  the  linear 
coefficient  values  obtained  when  Case  1 and  Case  2 were  combined  at 
zero  knot.  The  angular  rates  retained  relatively  unchanged  in  this 
situation.  However,  in  the  high  velocity  pooling,  both  the  linear 
and  angular  coefficient  values  were  reduced  from  their  corresponding 
individual  component  values  found  in  Tables  3 and  4.  This  difference 
can  be  attributed  to  the  differences  existing  between  the  two  sets  of 
input  parameter  base  values  used  in  the  study  (Case  1 versus  Case  2) 
and  specifically  to  the  large  difference  (20%)  in  the  seat/man  weight 
existing  between  the  two  cases  at  the  600  knot  initial  velocity  condition, 
whereas  in  the  0-knot  pooling  of  the  rankings,  a five  percent  difference 
between  the  two  corresponding  seat/man  weights  was  not  so  noticeable. 

The  between  case  pooling  of  the  individual  component  rankings  of  the  input 
parameters  indicates  that  the  rankings  of  the  input  parameters  can  be 
effected  by  differences  in  the  initial  values  of  the  parameters. 

In  several  instances  the  Spearman's  Rank  Correlation  Coefficients  in  Table 
8 were  slightly  higher  than  their  corresponding  component  values 
found  in  Tables  3 and  4.  Again,  this  was  due  to  the  numerical  averaging 
effect  inherent  in  the  ranking  procedure  used. 
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In  all  cases,  the  interpretation  of  the  Coefficients  of  Concordance, 
as  to  the  overall  degree  of  agreement  in  the  rankings  obtained,  was 
found  to  be  statistically  significant  at  the  95%  confidence  level. 

Tables  R9  and  RIO  present  the  actual  component  rankings  of  the  eight 
input  parameters  obtained  by  combining  the  Case  1 and  Case  2 corres- 
ponding individual  component  rankings.  Here  too,  the  rankings  formed 
by  pooling  are  statistically  less  significant  than  the  rankings  dis- 
played for  the  individual  components . (Tables  1 thru  4) . 


5.3.2  Between  Case  Linear  and  Angular  Component  Poolings.  Table  9 
presents  the  between  case  linear  and  angular  component  pooling  (x,  y 
and  z Case  1 combined  with  x,  y and  z Case  2:  p,  q and  r Case  1 com- 
bined with  p,  q and  r Case  2) . The  coefficient  values  shown  in  Table  9 
tend  to  be  lower  in  value  than  their  corresponding  individual  pooled 
component  values  found  in  Tables  7 and  8.  These  lower  values  can  be 
attributed  in  some  part  to  the  effects  of  data  averaging  but,  from  the 
prior  results  noted  for  the  between  and  within  case  component  poolings, 
there  is  reason  to  suspect  that  the  reductions  are  also  due  to  the  differ- 
ences existing  between  the  components  in  the  ranking  of  the  parameters. 
That  is,  the  individual  component  rankings  are  different  from  one 
another. 

Tables  Rll  and  R12  display  the  actual  rankings  of  the  input  parameters 
obtained  by  the  between  case  pooling  of  the  linear  trajectory  component 
rankings  (x,  y and  z combined)  and  angular  trajectory  component  rankings 
(p,  q and  r combined) . These  rankings  are  statistically  less  significant 
than  their  corresponding  individual  component  rankings  displayed  in 
Tables  R9  and  RIO. 


5.3.3  Summary  of  the  Results  Obtained  from  the  Ranking  Procedure  used 
in  the  Stuciy]  In  conclusion,  the  preceding  analysis  of  the  individual 
and  pooled  rankings  of  the  eight  input  parameters  by  the  six  trajectory 
components  indicates  the  following: 

1)  the  individual  trajectory  component  rankings  of  the  eight 
input  parameters  are  consistent  (in  agreement)  from  phase 
to  phase  of  the  trajectory  as  indicated  by  the 

high  values  of  the  Concordance  Coefficients  obtained  for 
each  of  the  above  rankings  formed  addition,  these 

coefficients  were  found  to  be  st*  _ly  significant 

at  the  95%  confidence. 

2)  as  the  component  rankings  art  pooled  together,  linear  component 
rankings  combined  and  angular  component  rankings  combined,  the 
degree  of  agreement  as  to  the  ranking  of  the  parameters 
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decreases.  This  is  confirmed  by  the  decrease  in  value  of 
the  Coefficients  of  Concordance  and  Spearman's  Rank 
Correlation  Coefficients.  This  implies  that  the  individual 
trajectory  component  rankings  can  be  different  from  one 
another. 

3)  the  further  decrease  in  the  Concordance  Coefficient  and 
Spearman's  Rank  Correlation  values  when  the  trajectory 
component  rankings  of  the  parameters  were  pooled  across 
velocities  (zero  velocity  combined  with  600  knots  velocity) 
indicates  that  the  ranking  of  the  parameters  obtained  at 
zero  velocity  was  different  than  the  ranking  of  the  para- 
meters at  600  knots  velocity. 

4)  the  decrease  in  the  values  obtained  for  the  Concordance 
Coefficients  and  Spearman's  Rank  Correlation  Coefficients 
when  the  individual  trajectory  component  rankings  of  Case 
1 were  combined  with  Case  2 indicates  two  physically 
different  sets  of  input  data  can  rank  the  variables  dif- 
ferently from  one  another. 

5)  the  weight  of  the  seat/man  was  identified  as  being  a sig- 
nificant factor,  in  addition  to  the  initial  velocity,  in 
determining  the  ranking  of  the  input  parameters. 

An  examination  of  the  plots  (Appendix  A)  made  of  Case  1 varied  parameter 
trajectories  versus  their  corresponding  base  case  trajectory  and  the 
tabulations  of  the  residual  values  (Appendix  B)  indicates  as  the  tra- 
jectory progresses  (from  catapult  ignition  to  1 second) , the  difference 
between  the  varied  parameter  trajectory  and  the  base  case  trajectory 
increases . 


5.3.4  Explanation  of  the  Coded  Rank  Tabulations  of  the  Input  Parameters. 
Tables  Dl  and  D2,  contained  in  Appendix  D,  present  two  numerically  coded 
tabulations.  Table  Dl  displays  a sunmary  of  the  individual  linear  com- 
ponent rankings  of  the  input  parameters  and  Table  D2  a sumrary  of  the 
angular  component  rankings  of  the  parameters  in  coded  form.  The  above 
coded  rankings  were  generated  using  Tables  Rl  and  R3  (Ref.  Appendix  C) 
to  obtain  the  linear  rankings  and  Tables  R2  and  R4  to  obtain  the  angu- 
lar rankings  of  the  parameters.  Tables  Rl  thru  R4  present  the  individual 
trajectory  component  rankings  of  the  parameters. 

The  format  of  the  coded  tabulation  being  an  8x8  matrix.  Each  of  the 
eight  columns  of  the  matrix  being  assigned,  and  titled  as  such,  to  one 
of  the  eight  input  parameters  selected  to  be  studied.  The  rows  of  the 
matrix  (numbered  1 to  8)  being  used  to  designate  the  relative  rank  of 
the  parameters. 
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The  numbers  1 thru  4 found,  in  the  individual  cells  of  the  two  matrices, 
designate  the  following: 

1.  the  number  1 implies  Case  1 at  0-knots  initial  velocity  pooled 
+ 20%  variation. 

2.  the  number  2 implies  Case  1 at  600-knots  initial  velocity  cooled 
+20%  variation. 

3.  the  number  3 implies  Case  2 at  0-knots  initial  velocity  pooled 
+ 20%  variation. 

4.  the  number  4 implies  Case  2 at  600-knots  initial  velocity  pooled 
+20%  variation. 

Replicates  of  the  above  numbers  contained  in  any  of  the  row/column 
designated  cells  of  the  two  matrices  indicating  the  number  of  times 
that  particular  parameter  (designated  by  the  column  heading)  was 
assigned  the  rank  indicated  by  the  row  number  of  the  cell  by  either 
the  x,  y or  z linear  trajectory  component  (Table  Dl)  or  the  p,  q or  r 
angular  trajectory  component  (Table  D2) . The  rank  of  the  specific 
parameter,  relative  to  the  other  seven  parameters,  being  determined 
by  the  amount  the  parameter,  when  varied  from  its  initial  base  case 
value,  caused  the  trajectory  component  value  to  deviate  from  its  base 
case  value. 

Tables  Dl  and  D2  present  a visual  summary  of  the  results  obtained  by 
varying  the  eight  input  parameters  individually  first  by  a plus  20%, 
then  by  a minus  20%  while  holding  the  other  parameters  constant  at  their 
respective  base  case  values.  An  inspection  of  the  two  tables  reveals 
the  following: 

1)  Six  of  the  eight  input  parameters  were  ranked  first 

(at  least  once)  by  one  of  the  six  trajectory  oorp  cnmts . 

The  other  two  weie  either  renked  second  or  third  one 
or  more  times.  Specifically,  the  moments  of  inertia 
were  ranked  second  once  by  one  of  the  angular  rates 
as  indicated  in  Table  D2.  The  Rocket  On  parameter 
was  ranked  third  three  times,  once  by  one  of  the  tra- 
jecotry  linear  coordinates  and  twice  by  one  of  the 
trajectory  angular  rates  as  indicated  in  Tabels  Dl 
and  D2  respectively. 

2)  the  rankings  of  the  parameters  by  the  linear  components 
(x,  y and  z)  and  the  angular  rate  components  (p,  q and 
r)  are  different  from  one  another. 

3)  the  rankings  of  the  parameters  are  effected  by  physical 
differences  in  the  ejection  systems  being  simulated. 

4)  the  rankings  of  the  parameters  are  effected  by  the  initial 
velocity  of  the  seat/man  combination. 
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5.4  PARAMETERS  VARIED  SIMULTANEOUSLY 

An  additional,  albeit  very  limited,  analysis  of  two  then  three  of  the 
input  parameters  being  varied  simultaneously  was  performed’.-  All  com- 
binations of  the  parameters  being  varied  by  a plus  or  minus  twenty 
percent  and  their  associated  ejection  trajectories  were  generated. 

The  catapult  and  rocket  time/thrust  tables,  base  case  values  for  the 
Case  1 at  600  knots,  were  first  varied  simultaneously  followed  by  the 
catapult  and  rocket  time/thrust  tables  plus  the  seat/man  center  of 
gravity  being  varied  simultaneously. 

As  would  be  expected,  certain  combinations  of  the  input  parameters 
being  varied  simultaneously  tended  to  minimize  the  difference  between 
the  varied  parameter  trajectory  obtained  and  its  corresponding  base 
case  trajectory,  while  other  combinations  tended  to  maximize  the 
difference . 

Appendix  A,  (figures  A-19  thru  A-24  show)  the  trajectories  obtained  by  the 
simultaneous  varying  of  the  above  parameters  along  with  their  corresponding 
base  case  trajectory. 


5.5  SPECIAL  CASES  SHOWING  THE  EFFECT  OF  VARYING  ALL  EIGHT  PARAMETERS 

SIMULTANEOUSLY 

Two  special  cases  (Case  1 at  0-knot  and  Case  1 at  600-knots  initial 
velocity)  were  devis<_  whereby  the  eight  input  parameters  were  varied 
by  small  amounts  simultaneously.  The  two  trajectories  obtained  (one 
for  the  0-knot,  the  other  for  the  600-knots  initial  velocity  condition) 
by  varying  the  eight  parameters  by  small  amounts  deviated  from 
their  corresponding  base  case  trajectories  almost  as  much  as  the 
maximum  deviation  obtained  by  varying  the  parameters  individually 
by  + 20%. 

Appendix  A,  (figures  A-25  and  A-26)  show  the  above  varied  parameter 
trajectories  along  with  their  respective  base  case  trajectory. 
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, ' RECOMMENDATIONS 

Because  of  the  near  infinite  number  of  possible  combinations  of  the 
input  parameter  values  attainable  by  varying  them  simultaneously  any- 
where from  zero  to  a plus  or  minus  twenty  percent,  an  exhaustive 
case  by  case  study  of  the  parameter  interaction  terms  would  be 
impractical.  If  a more  indepth  analysis  of  the  simulation  model's 
response  to  varying  specific  or  all  of  the  input  parameters  simul- 
taneously by  the  same  or  different  amounts  is  required,  a stochastic 
approach  might  be  considered.  The  fact  that  all  of  the  input  parameters' 
base  case  values  have  seme  degree  of  uncertainty  associated  with  them 
supports  this  type  of  approach  being  taken. 


6.1  OUTLINE  OF  THE  STOCHASTIC  APPROACH 

The  stochastic  approach  was  not  used  for  this  effort  due  to  very  high 
projected  computer  cost,  however,  it  would  be  possible  to  allow  the  user 
of  the  model(s)  to  enter  the  parameters  to  sampling  distributions.  These 
distributions  would  then  be  used  to  obtain  the  model's  input  parameter 
values.  The  actual  values  of  the  input  parameters  being  randomly  gen- 
erated from  their  respective  sampling  distribution.  An  ejection  trajectory 
would  then  be  generated,  via  the  model,  based  upon  these  randomly  gen- 
erated values  of  the  input  parameters.  The  random  input  parameter 
values  and  ensuing  trajectory  generation  would  be  repeated  n- times.  After 
n-iterations,  an  envelope  of  trajectories  surrounding  the  base  case 
trajectory  would  be  generated.  Once  the  proper  sample  size  (n)  is  ascer- 
tained, the  boundaries  (or  near  boundaries)  of  the  envelope  should  be 
among  the  (n)  trajectories  generated. 

Thus,  based  upon  the  user's  best  estimate  as  to  the  degree  of  variability 
existing  around  each  input  parameter,  the  maximum  expected  difference 
between  the  base  case  trajectory  (generated  using  the  most  likely  values 
of  the  input  parameters)  and  the  n-randomly  generated  trajectories 
should  be  ascertained. 

Figures  1 to  4 display  the  envelope  surrounding  each  of  the  four  base 
case  trajectories  (Case  1 and  2 at  zero  and  500  knots  respectively). 

The  envelopes  reflect  the  bounds  of  the  trajectories  generated  by  varying 
the  parameters  individually  by  a plus  or  minus  20%. 
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SECTION  7 

DOCUMENTATION  OF  THE  LATEST  REVISIONS  TO  THE  ICARUS  AND  DEDALUS  SEAT 

EJECTION  MODELS 


The  following  documents  the  work  completed  to  date  by  Computer  Sciences 
Corporation  in  modifying  and  extending  the  capabilities  of  two  seat 
ejection  mathematical  simulation  models,  namely: 

1)  ICARUS,  a 3-dimensional  model  which  completely  describes  the 
ejection  trajectory  of  the  seat/man  combination  by  computing 
at  each  time  step  of  the  trajectory . the  3 spatial  coordinates 

(X,  Y and  Z)  and  the  3 angles  of  orientation  (roll,  pitch  and  yaw)  of 
the  seat/man  combination  relative  to  an  earth  fixed  coordinate  system. 

2)  DEDALUS,  a 2-dimensional  model,  which  describes  the  ejection 
trajectory  of  the  seat/man  combination  by  computing  at  each 
time  step  of  the  trajectory  the  two  linear  coordinates  (X 
coordinate-horizontal,  Z coordinate-vertical)  and  the  pitch 
angle  in  the  plane  of  X and  Z coordinates. 


7.1  MODIFICATIONS  AND  EXTENSIONS  TO  ICARUS'  AERODYNAMIC  COEFFICIENT 
TABLES 

The  ICARUS  program  is  designed  to  use  seat/man  combination  aerodynamic 
coefficients  contained  in  a random  file.  The  coefficients  previously 
used  were  obtained  from  Aeromechanical  Properties  of  Ejection  Seat  Escape 
Systems.  These  coefficients  were  used  in  their  raw  form  as  given  by  the 
above  report.  A more  thorough  analysis  of  these  data  has  been  performed 
resulting  in  a more  consistent  and  expanded  set  of  seat/man  combination 
aerodynamic  coefficients  for  use  by  the  ICARUS  program. 


7.1.1  Smoothing  of  the  Aerodynamic  Angles  of  Attack.  Due  to  the  limi- 
tations of  the  sting  type  mounting  system  used  in  the  wind  tunnel  test, 
it  was  only  possible  to  test  a range  of  angles  of  attack  of  0°  to  120° . 
However,  by  changing  the  attachment  point  of  the  sting  to  the  seat/man 
model,  the  angle  of  attack  range  of  120°  to  240°  could  be  obtained. 

In  the  same  manner  the  angle  of  attack  range  of  240°  to  360°  could 
also  be  tested.  When  these  data  are  plotted  versus  angle  of  attack 
the  influence  of  the  new  sting  location  relative  to  the  seat/man  model 
shows  up  as  a jump  in  coefficient  magnitude  at  angles  of  attack  of  0°, 
120°,  and  240°.  Since  the  true  coefficient  values  should  be  between 
these  two  levels  of  magnitude  a median  fairing  line  was  drawn  between 
the  coefficients  at  these  angles  of  attack.  Also,  a best  approximation 
fairing  line  was  drawn  through  all  the  raw  data  points  for  the  full  angle 
of  attack  range.  This  analysis  results  in  a smoothed  set  of  aerodynamic 


1.  White,  B.J.,  "Aeromechanical  Properties  of  Ejection  Seat  Escape 
Systems",  AFFDL-TR-74-57 , April  1974 
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coefficients  and  these  curves  were  then  read  at  the  desired  angles  of 
attack  to  obtain  the  aerodynamic  coefficients  versus  angle  of  attack. 
These  tabulated  points  were  then  plotted  versus  Mach  number  and  angle 
of  sideslip.  These  plots  were  then  read  at  the  desired  Mach  numbers  and 
angles  of  sideslip  to  yield  the  complete  set  of  aerodynamic  coefficients 
to  be  used  by  the  ICARUS  program. 


7.1.2  Extension  of  ICARUS  Aerodynamic  Sideslip  Coefficient  Tables  for  Sideslip. 

The  previous  version  of  the  ICARUS  program  was  not  designed  to  handle 

aerodynamic  coefficients  for  a full  range  of  sideslip  angles.  Since 

the  raw  data  frcm  the  wind  tunnel  test  only  included  sideslip  angles 

of  0°  to  45°  an  estimate  was  made  for  the  aerodynamic  coefficients 

for  sideslip  angles  of  45°  to  90°.  If  sideslip  angles  greater  than 

90°  were  generated  by  the  ICARUS  program  the  run  was  terminated. 

The  aerodynamic  coefficients  from  the  wind  tunnel  test  are  in  the 
Seat/Man  Coordinate  System  which  are  fixed  to  the  seat/man  combina- 
tion. Therefore,  the  aerodynamic  coefficient  at  any  angle  of  attack  (a) 
and  sideslip  angle  ( ip  ) is  duplicated  at  the  following  orientation : 

oc'=  a + 180° 

ip'=  180°  - ip 

Using  this  relationship,  a plot  of  coefficients  versus  sideslip  angles  of  0°  to 
45°  and  135°  to  130°  can  be  made.  This  plot  can  then  be  interpolated  to 
obtain  coefficients  in  the  45°  to  135°  sideslip  angle  range. 

Due  to  the  limitations  of  its  design  the  ICARUS  program  will  never 
generate  a sideslip  angle  greater  than  90° . If  the  ICARUS  simulation 
run  physically  tries  to  rotate  past  a sideslip  angle  of  90°  the  angle 
of  attach  and  sideslip  angle  are  calculated  so  as  to  have  the  same 
coefficient  as  if  it  had  a sideslip  angle  greater  than  90° . Therefore, 
it  was  only  necessary  to  expand  the  aerodynamic  coefficient  tables  to 
include  sideslip  angles  of  0°  to  90° . This  was  done  using  data  from 
the  previously  described  plots  of  coefficients  versus  sideslip  angles. 

The  subroutine  SEATMAN  was  changed  to  remove  the  restrictions  on  side- 
slip angle  since  there  are  none  with  the  use  of  the  modified  aerodyna- 
mic coefficient  tables.  It  was  also  necessary  to  modify  subroutine 
ACT  in  order  to  handle  the  increased  size  of  these  expanded  aerodyna- 
mic coefficient  tables. 


7.2  MODIFICATIONS  AND  EXTENSIONS  TO  ICARUS'  ATMOSPHERIC  DENSITY, 

SPEED  OF  SOUND,  AND  MACH  NUMBER  CALCULATIONS. 

The  ICARUS  program  was  modified  so  that  it  would  calculate  the  atmos- 
pheric density,  speed  of  sound,  and  object  Mach  number  for  a specified 
pressure  and  temperature  rather  than  the  previously  assumed  standard 
day  conditions. 
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The  atmospheric  conditions  are  defined  by  the  ambient  pressure  and 
temperature  existing  at  the  time  and  point  of  the  ejection.  The  sub- 
routine INPOT  was  modified  to  accept  the  pressure  (PMBS)  in  milibars 
and  the  temperature  (TFAR)  in  °F  as  items  13  and  14  of  section  3 of 
the  input  data  file.  Subroutine  INIT  was  expanded  to  calculate  con- 
stant parameter  values  based  on  the  given  temperature  and  pressure. 
Subroutine  SEA.TMAN  was  changed  to  calculate  the  atmospheric  density, 
speed  of  sound,  and  object  Mach  number  for  each  iteration  of  the 
ejection  trajectory. 


The  constant  values  calculated  by  subroutine  INIT  are  as  follows: 


Given:  Pamb  = Ambient  pressure  in  milibars 

Tamb  = Ambient  temperature  in  degrees  Fahrenheit  (°F) 

Calculate:  T'q^  = Ambient  temperature  in  degrees  Rankine  (°R) 

= T^  + 459.688 


Hp  = Pressure  altitude  in  feet  for  Pamb 

/ Pamb  \ 0.19 
_ 1 \1013.25/ 

0.00006875 


When  the  above  Hp  >36089  feet  it  must  be 
recalculated  by: 


Hp  = -In 


Pamb  \ \ 

05#  / 20787  + 36089 


Tstd  = Standard  day  temperature  in  °R  at  Hp 
= 518.688  - (0.003566  Hp)  when  Hp  £36089 
= 390.0  when  Hp  >36089 

AT  = Temperature  deviation  frcm  standard  in  °R  at  Hp 
= T'amb  ~ Tstd 


For  each  iteration  during  the  ejection  trajectory  subroutine  SEATMAN 
calculates  the  corresponding  atmospheric  density,  speed  of  sound,  and 
object  Mach  number  as  follows: 

Given:  V - Current  object  velocity  in  feet/second 

Calculate:  AH  = Change  in  trajectory  height  in  feet 
= H - Hlast 

Hp  = Current  pressure  altitude  in  feet 
= Hp  + AH 
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Tstd  = Standard  day  temperature  in  °R  at  Hp  assuming  a 
standard  day  lapse  rate 
= 518.688  - (0.003566  Hp)  when  Hp  <36089 
= 390.0  when  Hp>  36089 
T'amb  = Ambient  temperature  in  °R  at  Hp 
= Tstd  + AT 


pstd  = Standard  day  pressure  in  milibars  at  Tst^ 
Tstd  \ 5-256  when  Hp£  36089 


= 1013.25  

518.688 

_ 1013.25  x 0.2234 


/ Tstd  \ 
\ 518.688  / 


2.71828^ 
where  X 


when  Hp  > 36089 


_ Hd  -36089 
20787 

= Atmospheric  density  in  slugs/feet^  at  Hp  and  Tg^ 

=0.0023769  / pstd  \ / 518.688  \ 

\ 1013.25  ) [ T^  j 


Vs  = Speed  of  sound  in  feet/second  at  Hp  and 
= 49.0212 

M = Mach  number  of  the  object  at  Hp,  Tamb>  and  V 
_ V 


7- 3 MODIFICATIONS  AND  EXTENSIONS  TO  ICARUS'  SEAT/MAN  INITIAL  PARTIAL 
EXPOSURE  CALCLTATIONS. 

7.3.1  Definition  of  the  Partial  Exposure  Input  Paramters.  The  ICARUS 
program  has  been  modified  to  use  additional  input  parameters  in  an  effort 
to  more  accurately  define  the  aerodynamic  forces  and  moments  on  the  seat/ 
man  combination  while  it  is  moving  out  of  the  cockpit  during  the  initial 
phase  of  the  ejection  sequence.  Parameters  were  added  to  account  for 
the  shielding  effect  of  the  windscreen  or  any  piece  of  equipment,  the 
exposed  area  of  the  seat/man  combination  due  to  crew  position  within 
the  cockpit,  and  the  varying  aerodynamic  forces  and  moments  acting  on 
the  seat/man  combination  during  the  ejection  sequence  until  the  cockpit 
has  beem  cleared. 

Subroutine  INPUT  was  modified  to  accept  the  following  parameters  as  items 
35  through  40  of  section  7 of  the  input  data  file. 


REPORT  IIO.  MAD 07 7 100-40 


CREW  - Crew  member  being  ejected  ) 

0 = Pilot 

1 = MCO  I 

TRAVEL  - The  height  (bottom  of  boots  to  top  of  seat)  of  the  seat/man 
combination  in  feet. 

SCRENO  - The  distance  in  feet  that  the  pilot  has  to  travel  before 

being  exposed  to  the  airstream.  ] 

SCREN1  - The  distance  in  feet  that  the  MCO  has  to  travel  before 

being  exposed  to  the  airstream.  j 

SUN  - The  exposed  area  in  square  feet  of  the  seat/man  combi- 
nation installed  in  the  cockpit. 

CPCG  - This  parameter  defines  the  point  at  which  the  transition  ' 

is  made  from  estimated  aerodynamic  moments  to  fully  1 

exposed  aerodynamic  moments.  It  is  the  ratio  of  partially 
exposed  area  to  fully  exposed  area  at  which  the  aerodyna- 
mic pitching  moment  will  be  0.  j 

Also,  subroutine  INPOT  was  modified  to  accept  the  following  parameters 
as  items  74  through  76  of  section  2 of  the  input  data  file. 

XCP,  YCP,  ZCP  - The  X,  Y,  and  Z coordinates  in  inches  of  the 

estimated  center  of  pressure  applicable  when  I 

the  seat/man  combination  first  enters  the  1 

airstream. 

! 

Refer  to  Figure  5 for  a description  of  these  parameters.  j 

7.3.2  Explanation  of  the  Partial  Exposure  Logic  and  Equations  Imple-  i 

men ted  in  ICARUS.  As  the  ejection  sequence  is  executed  the  seat/man 

combination  will  progress  from  its  initially  shielded  or  partially  I 

exposed  position  until  it  is  fully  immersed  in  the  airstream.  Once 

this  occurs,  the  fully  exposed  aerodynamic  coefficient  and  areas  may  i 

be  used  to  calculate  the  aerodynamic  forces  and  moments  acting  on  the 
seat/man  combination.  However,  prior  to  this  event  it  will  be  neces- 
sary to  estimate  the  aerodynamic  forces  and  moments  acting  on  the  seat/ 
man  combination  while  it  is  influenced  by  the  cockpit  configuration, 

since  the  fully  exposed  aerodynamic  coefficients  and  areas  are  not  J 

applicable  to  this  situation. 

' I 

At  the  time  of  ejection  initiation  the  seat/man  combination  will  be 
shielded  from  the  airstream  or  will  be  exposed  to  some  extent  depend- 
ing on  the  crew  position  and  cockpit  configuration.  For  the  shielded 
case  (SCRENO  or  SCRENl)  the  seat/man  combination  must  move  a finite 
distance  before  it  enters  the  airstream  and  aerodynamic  forces 
and  nonents  are  developed.  However,  the  partially  exposed  (SUN) 
seat/man  combination  already  has  aerodynamic  forces  and  moments  acting 
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on  it  at  the  time  of  ejection  initiation.  Therefore,  for  either  case 
the  catapult  tubes  must  extend  some  fixed  distance  for  the  seat/man 
combination  to  be  fully  immersed  in  the  airstream.  This  distance 
(CLEAR)  is  determined  in  subroutine  INIT  by: 

CLEAR  = SCREN  + TRAVEL  - 

WIDE 

where,  SCREN  = SCRENO  or  SCRENl 

WIDE  = ss 

TRAVEL 

where,  SS  is  the  fully  exposed  area 

Until  the  catapult  tubes  have  extended  to  equal  CLEAR  the  partially 
exposed  area  will  be  estimated  by: 

SEX  = SUN  + ( (TUBEXl  - SCREN)  WIDE) 

where,  TUBEXl  is  the  extension  of  the  tubes 

SCREN  = SCRENO  or  SCRENl 

The  above  estimated  areas  will  be  calculated  by  subroutine  AEKOEM  in 
order  to  estimate  the  aerodynamic  forces  acting  on  the  partially  exposed 
seat/man  combination  until  full  exposure  has  been  achieved.  The  fully 
exposed  aerodynamic  force  coefficients  will  be  used  regardless  of  ex- 
posure. 

The  aerodynamic  moments  acting  on  the  seat/man  combination  while  partially 
exposed  will  be  approximated  by  using  the  above  estimated  aerodynamic 
forces  acting  over  an  estimated  moment  arm.  The  initial  estimated  moment 
arms  are  calculated  by  subroutine  INIT  as  the  difference  between  the 
estimated  center  of  pressure  coordinates  (XCP,  YCP,  ZCP)  and  the  seat/ 
man  center  of  gravity  coordinates  (XCS,  YCS,  ZCS) . As  the  seat/man 
combination  moves  out  into  the  airstream  the  Z moment  arm  will  be 
assumed  to  decrease  as  follows: 

ZARM  = ZARM  - | ZARM 

where,  SX  = SEX 

35^ 

From  the  above  equation,  the  Z moment  arm  will  be  reduced  to  0 when  SX  = 
CPCG.  Beyond  this  point  the  fully  exposed  aerodynamic  moment  coefficient 
will  be  used. 

Subroutine  AEROFM  checks  to  see  if  the  tube  extension  is  greater  than 
SCREN  (SCRENO  or  SCRENl) . If  not,  then  the  exposed  area  is  set  to  0 
and  therefore  aerodynamic  forces  and  moments  of  0 will  be  calculated. 
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When  the  tube  extension  is  greater  than  SCREN  a check  is  made  to  see  if 
the  tube  extension  is  greater  than  CLEAR.  If  so,  the  fully  exposed  area 
is  used.  Otherwise,  the  estimated  exposed  area  is  calculated  and  used. 

The  following  input  parameter  values  are  illegal  conditions  and  will 
result  in  an  error  message  and  termination  of  the  run: 

(SCRENO  or  SCRENl)  >0  and  SUN  >0 

(YCP  - YCS)  >TRAVEL 

CLEAR  > (TUBLTHl  + TUBLTH2) 

(SCRENO  or  SCRENl)  > TRAVEL 

SUN  >SS 

CPCG  >1.0 


7.4  REVISIONS  TO  BOTH  DEDALUS ' AND  ICARUS'  SIMULATION  OF  THE  MARTIN-BAKER 
CATAPULT  TUBE  BENDING  AND  SEAT  BOTTOM  FORCES. 

The  following  documents  the  latest  revisions  made  to  both  the  DEDALUS  and 
ICARUS  seat  ejection  models  with  respect  to  the  simulation  of  the  Martin- 
Baker  Catapult. 

7.4.1  Replacement  of  the  Martin-Baker  Tube  Bending  Equations  with  a 
Spring  Equation.  The  tube  bending  equation,  originally  entered  into 
both  the  DEDALUS  and  ICARUS  seat  ejection  models  to  simulate  the  forces 
and  moments  acting  on  the  seat/man  combination  as  the  catapult  tubes 
extend,  have  been  replaced  with  a spring  equation  of  the  form 

Ftube  = C1  " (r/x)P  } {kAD  + CAD} 

where,  AD  is  the  resultant  distance  that  the  tube  has  bent  forward 
or  backward  (+X  - direction)  and  moved  sideways  (yawed,  + Z 
direction)  in  the  catapult  coordinate  system  relative  to 
its  initial  static  position. 

R is  the  length  the  tube(s)  have  extended  since  catapult 
ignition. 

X is  the  total  length  of  the  fully  extended  catapult, 
p is  an  empirically  fitted  oonstant. 
k is  the  stiffness  constant  (pounds/feet) 
c is  the  danping  coefficient  (pounds-second/feet) 


-33- 


REPORT  NO.  NAD C-77 100-40 


Ftube  is  the  resultant  force  (lbs.)  acting  on  the  seat/man 
con±)ination  due  to  the  tube  bending. 

A friction  force  due  to  the  tube  bending  is  also  computed  as  well  as  the 
components  of  the  Ftube  force.  The  tube  bending  force  is  computed  in 
the  catapult  coordinate  system  and  transformed  to  the  seat/man  coordi- 
nate system.  The  calculation  of  this  force  is  ended  once  tube  separation 
occurs. 


7.4.2  Insertion  of  a Seat  Bottom  Restoring  Force.  A restoring  force  is 
computed  at  the  bottom  of  the  seat  which  prevents  the  seat  from  moving 
through  the  tube.  If  X is  the  distance  that  the  bottom  of  the  seat  has 
traversed  through  the  tube  in  the  plane  perpendicular  to  the  tube  then 
the  magnitude  of  this  force,  F,  is  given  by 

F = -KX 

where  K is  the  stiffness  constant.  The  calculation  of  this  force  is  com- 
pleted once  the  bottom  slippers  have  cleared  the  top  of  the  rail. 


7.5  MODIFICATIONS  AND  EXTENSIONS  TO  DEDALL'S'  ATMOSPHERIC  DENSITY,  SPEED 
OF  SOUND  AND  MACH  NUMBER  CALCULATIONS. 


The  DEDALUS  program  has  been  modified  to  include  the  effect  of  the  local 
pressure  and  temperature  on  the  attnosphereic  density  and  speed  of  sound 
used  in  the  aerodynamic  force  and  moment  calculations. 


Subroutine  INPUT  was  altered  to  accept  the  following  atmospheric  para- 
meters from  the  input  data  file: 

- Atmospheric  temperature  in  degrees  Fahrenheit  (°F) 

PMBS  - Atmospheric  pressure  in  milibars 


Subroutine  INPUT  was  changed  so  that  the  atmospheric  density  and  speed 
of  sound  will  be  calculated  based  on  the  input  pressure  and  temperature. 
The  atmospheric  density  and  speed  of  sound  remain  constant  throughout 
the  trajectory  and  are  not  corrected  for  changes  in  altitude.  The 
equations  are: 


I’amb  = Ambient  temperature  in  degrees  Rankine  (°R) 
= + 459.688 


V 


s 


P 


= Local  speed  of  sound  in  feet/second 
= 47.02 


= Atmospheric  density  in  slugs/feet 


3 


= 0.002379 


/ PMBS  \ 
^1013.25 j 


I 518.688  \ 
\ Tamb  / 
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7.6  REPLACEMENT  OF  DEDALUS  AERODYNAMIC  COEFFICIENT  CURVE  FITTED  EQUA- 
TIONS WITH  AERODYNAMIC  COEFFICIENT  TABLES. 

The  method  of  providing  the  seat/man  combination  aerodynamic  coefficients 
for  the  DEDALUS  program  has  been  modified.  Previously,  subroutine  AERSEAT 
contained  curve  fitted  equations  based  on  the  raw  test  data  from  Aerome- 
chanical  Properties  of  Ejection  Seat  Escape  Systems.  These  curve  fitted 
equations  were  then  used  to  calculate  the  coefficients  for  a given  Mach 
number  and  angle  of  attack.  Due  to  the  errors  inherent  in  trying  to 
curve  fit  these  aerodynamic  coefficient  data  it  was  decided  to  revert  to  a 
table  presentation  of  the  aerodynamic  coefficients.  The  curve  fitted 
equations  of  subroutine  AERSEAT  were  replaced  with  DATA  statements  con- 
taining the  coefficients  as  a function  of  angles  of  attack  and  tiach 
numbers.  These  tables  are  interpolated  to  obtain  the  coefficients 
for  the  desired  Mach  number  and  angle  of  attack. 


7.7  IMPLEMENTATION  OF  A SUMMARY  OUTPUT  REPORT  INTO  DEDALUS. 

The  DEDALUS  program  was  modified  to  provide  the  printing  of  a summary 
report . The  existing  trajectory  report,  aerodynamic  report  and 
catapult  reports  can  now  be  selected  as  options. 

The  sunmary  report  will  always  be  output  and  consists  of  events  only. 
Each  time  an  event  occurs  the  appropriate  message  is  printed  along  with 
time,  horizontal  distance,  vertical  distance,  pitcn  angle,  Mach  number, 
horizontal  velocity,  and  vertical  velocity.  The  optional  reports  can 
be  printed  by  setting  the  following  input  data  file  flags  to  1. 

INTRAJ  - Print  trajectory  report 
IAERO  - Print  aerodynamic  report 
IOPT  - Print  catapult  reports 

Setting  any  of  these  flags  to  0 will  surpress  that  report. 


7.8  IMPLEMENTATION  OF  A SUMMARY  OUTPUT  REPORT  INTO  ICARUS 

The  ICARUS  program  was  modified  to  provide  the  user  of  the  model  with  the 
option  of  printing  a summary  output  report  in  addition  to  or  in  lieu  of 
i the  existing  detailed  trajectory  output  reports.  The  report  displays 

' time,  the  components  of  the  seat/nan  acceleration,  velocity,  linear 

i i position  and  angular  velocity  and  the  rotation  angles  used  to  rotate 

from  the  seat/man  coordinate  system  to  the  earth  fixed  coordinate  system. 
The  above  variables  are  displayed  in  truncated  form  with  1 digit  after 
the  decimal  point  being  retained. 

Location  25  of  the  IRSW  vector  contained  in  the  CARDIN  file  is  used  to 
‘ control  the  generation  or  to  inhibit  the  generation  of  this  report. 

Setting  IRSW  (25)  to: 
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0 - will  inhibit  the  generation  of  the  summary  report  while  having  no  effect 

on  the  other  reports  requested  by  the  user  of  the  model. 

1 - will  cause  the  model  to  generate  and  output  the  summary  report 

and  to  inhibit  the  generation  of  the  detailed  reports  if  any  have 
been  requested. 

2 - will  cause  the  summary  report  to  be  generated  and  outputted  in  addition 

to  the  other  reports  requested  by  the  user  of  the  model. 
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♦Indicates  the  significance  of  the  coefficient  of  concordance. 
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Multiple  entries  in  row/ column  indicates  same  rank  assigned  to  parameters  by  ranking 
procedure  used. 
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Multiple  entries  in  a row/column  indicate  the  sane  rank  was  assigned  to  parameters 
ranking  procedure  used. 
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CODED  RANK  TABULATIONS  OF  THE  INPUT  PARAMETERS 
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MATHEMATICAL  PRESENTATION  OF  THE  RESIDUAL  VALUES 

The  following  equations  were  used  to  calculate  the  average  value  of 
the  residuals  at  the  95%  confidence  level. 

1.  |res|  = |b|  - | vp | 

2.  RES^  j,k  = ®i,j,k  — VEi,j  k 

2 Ni  2 

3.  ^ k = s i k ~ RES.  , ) /N-; 

i=l  i/3»K  3,ky  3 

Where: 

| B | is  a 3-dimensional  matrix  containing  the  values  of  the 

six  seat  ejection  trajectory  components  (the  3 linear 
coordinates  x,  y and  z and  the  3 angular  rates  p,  q and 
r) . These  trajectory  components  were  obtained  by  simulat- 
ing the  ejection  of  the  seat/nan  via  ICARUS  with  the  input 
parameters  set  at  their  base  case  values. 

| VP | is  a 3-dimensional  matrix  containing  the  values  of  the  six 
seat  ejection  trajectory  components  (the  3-linear  coordi- 
nates x,  y and  z and  the  3-angular  rates  p,  q and  r) . 

These  trajectory  components  were  obtained  by  Simula ting 
the  ejection  of  the  seat/man  via  ICARUS  with  one  or  more 
of  the  eight  select  input  parameters  varied  from  their 
base  case  values,  while  the  other  input  parameters  were 
held  constant  at  their  base  case  values . 

| RES | is  the  3-dimensional  matrix  of  residual  values.  This 

matrix  was  obtained  by  subtracting  each  of  the  6-trajectory 
component  profiles  generated  for  the  varied  parameter 
case  from  their  corresponding  base  case  profiles. 

RESj^k  is  the  2-dimensional  matrix  containing  the  average  values 
' of  the  residuals  over  selected  intervals  of  the  seat/man 
ejection  trajectory. 

2 

Vj(k  is  the  2-dimensional  matrix  containing  the  variances 
computed  about  the  average  residual  values. 

i = 1,  2,  . . . , Nj ; where  = total  number  of  trajectory 
elements  in  the  jth  interval. 

j = 1 implies  the  interval  from  catapult  ignition  to  catapult 
separation. 

= 2 implies  the  interval  frcm  rocket  ignition  to  rocket 
burnout. 

= 3 implies  the  interval  from  rocket  burnout  to  drogue 
chute  projection. 
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= 4 implies  the  interval  from  drogue  chute  projection  to 
one  second  of  the  trajectory  simulation. 

k = 1,  2 and  3 implies  the  3-linear  coordinates  (x,  y and  z)  of 
the  seat/man  center  of  gravity  respectively  expressed  in  an 
earth  fixed  coordinate  system  in  units  of  feet. 

=4,5  and  6 implies  3-angular  rates  (p,  q and  r)  of  the 
seat/man  ccmbinates  with  origin  at  the  seat/man  center 
of  gravity,  expressed  in  radians/seconds. 

4.  Rj,k  = ^^jjk  + Cj  k 0j,k 

Where  0 < C j f ^ < 5 

Equation  4 was  used  to  obtain  the  2-dimensional  matrix 
of  residual  values  at  the  95%  confidence  level.  That  is, 
95%  of  the  residual  values  of  the  kth  trajectory  component 
contained  in  the  jth  trajectory  interval  where  < to  this 
value. 

The  above  Rj  ^ values  were  then  used  to  rank  the  eight 
input  parameters  selected  to  be  varied  in  the  study. 

. 
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APPENDIX  F 

MATHEMATICAL  EQUATIONS  USED  TO  COMPUTE  THE  COEFFICIENTS  OF 
CONCORDANCE  AND  SPEARMAN'S  RANK  CORRELATION  COEFFICIENTS 


■ X- 


Where:  F3\  = Rank  Total,  sum  of  the  individual  phase  ranking  values 

assigned  to  the  i^1  input  parameter. 

EV  = Expected  value  of  the  parameters 
^max  = Max™um  possible  sum  of  squares 

m = number  of  phases  or  intervals  used  to  rank  the  parameters 
n = number  of  items  to  be  ranked  (8) 

W = Coefficient  of  Concordance 
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F-value  - F-value  associated  with  the  Coefficient  of  Concordance 

Ravo  = average  value  of  Spearman ' s Rank  Correlation 
y Coefficient.  .. 

1=1,  2,  8 (Number  of  parameters  to  be  ranked) . 


F2 


* 


. 


DISTRIBUTION  LIST 

REPORT  NO.  NADC-77100-40 
AIRTASK  NO.  A5365360/001D/8W0625-0000 
Work  Unit  No.  A5313-34B 

No.  of  Copies 

NAVAIRSYSCOM,  AIR-950D  6 

(2  for  retention) 

(1  for  AIR-531) 

(3  for  AIR-53126 

DDC 12 

McDonnell  Douglas,  St.  Louis  (Charleville)  2 

Martin-Baker  Aircraft  Co.,  CTD  2 

Naval  Air  Test  Center  (SY-711) 2 

National  Parachute  Test  Range  2 


